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Abstract: We study the topological structure of the SU(2) vacuum at zero 
temperature: topological susceptibility, size, shape and distance distributions 
of the instantons. We use a cooling algorithm based on an improved action 
with scale invariant instanton solutions. This algorithm needs no monitoring 
or calibration, has an inherent cut off for dislocations and leaves unchanged 
instantons at physical scales. The physical relevance of our results is checked 
by studying the scaling and hnite volume dependence. We obtain a suscep¬ 
tibility of (200(15)MeV)'^. The instanton size distribution is peaked around 
0.43fm, and the distance distribution indicates a homogeneous, random spa¬ 
tial structure. 


1 Introduction 

Much work has been devoted to the study of topology within the framework of lattice 
gauge theories. Among other topics the measurement of the topological susceptibility and 
investigations of the relevance of instantons to the QGD vacuum have attracted a lot of 
interest. 

The topological susceptibility is a physical quantity which measures the fluctuations 
of the topological charge of the vacuum. It is related via the Witten-Veneziano formula 
to a phenomenological parameter (the rj' — rj mass difference) and attributed to the 
presence of instantons in the QGD vacuum. It can be measured on the lattice where it is 
expected to scale according to its dimension as the continuum limit is approached. 

Instantons are classical solutions of the (euclidean) equations of motion, that is, they 
represent minima of the euclidean action [^, They have been conjectured to be relevant 
for various effects like the hadronic spectrum or the chiral phase transition and used as 
the basic degrees of freedom of several effective models. Such models assume certain 
characteristics of the instanton ensemble concerning for instance the density correlations 
(crystal, liquid or gas), the size distribution and other features |^. Tests of the assumptions 
on which these models are based, or direct tests of the physical effects of instantons can be 
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provided by the study of the physical, non-trivial topological excitations present in lattice 
conhgurations. 

In connection with the measurement of the topological susceptibility on the lattice most 
approaches have concentrated on devising appropriate topological charge operators. The 
operators commonly used belong to three categories: 

• Geometric, based on a geometric construction of the topological charge 1,0- 

• Analytic, based on the held-theoretical definition of the charge 

• Fermionic, based on the relation between the topological charge and the zero modes 
of the Dirac operator ill, 00. 


The latter, although promising, has not been pursued too far yet due to the problems 
involved in taking the zero mass limit of the Dirac operator. There has also been a recent 
proposal for simulating chiral fermions on the lattice (the overlap formalism) [0 which 
simultaneously provides a method to measure the topological charge. It has been applied 
to measure the topological susceptibility and preliminary results indicate good agreement 
with the Witten-Veneziano formula |T| . 

The other two approaches face several difficulties related to lattice artifacts due in 
particular to the bad behavior of the Wilson action regarding scale-invariance. Instead 
of assigning to an instanton an action independent of its size, the Wilson action gives 
an action which monotonically decreases with the instanton size. This leads, in a Monte 
Carlo simulation, to an unphysical abundance of small instantons at the scale of the cut¬ 
off (“dislocations”). These dislocations are not physical but still they contribute to the 
geometrical charge which hence tends to overestimate the topological susceptibility 0, |6| . 
According to an argument by Pugh and Teper [W| their contribution will even diverge in 
the continuum limit. On the other hand naive field-theoretical definitions of the topological 
charge have strong renormalization effects induced by the lattice regularization which tend 
to decrease the topological susceptibility and again spoil its scaling behavior [0 . 

These problems can be handled at two levels: in the Monte Carlo simulation, by im¬ 
proving the scale invariance of the action used in order to decrease the contribntion from 
dislocations; or in the measurement, by designing a method to disentangle dislocations 
from physical instantons and to decrease renormalization effects. 

Both approaches have been pursued. Cooling (an iterative minimization of the lattice 
action) |T|| has been used to handle the problem at the level of measurement. Ideally 
it eliminates rough topological fluctuations (dislocations) while keeping large instantons 
unchanged and decreases renormalization effects by smoothing out the ultraviolet noise. 
For the configurations remaining after cooling good agreement is found between field- 
theoretical and geometric determinations of the topological charge. However cooling with 
the Wilson action indnces more changes than just smoothing out the UV noise and elimi¬ 
nating dislocations. The decrease of the Wilson action with the instanton size means that 
this action has actually no stable instantons: under cooling with the Wilson action even 
wide (physically relevant) instantons shrink and decay. Hence cooling algorithms based on 
such action cannot provide a reliable plateau for the topological charge as a function of 
cooling time. Moreover these algorithms do not preserve the physical scale structure of the 
original configurations. Because of this one usually attempts to calibrate or engineer the 
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cooling procedure based on the Wilson action so as to obtain some meta-stability window 
20|, after the noise has been reduced but before instantons have shrunk to zero. This, 


however, makes cooling more an art than a method and also does not unambiguously pre¬ 
serve information concerning, for instance, the physical size of the instantons present in 
the original configuration. 

There have been quite successful attempts to combine cooling with other methods, like 
“heating” |18], to compute the lattice renormalization factor of the topological suscepti¬ 


bility, or more recently to define “smeared” topological charge operators [^. They were 
shown to improve susceptibility data but no results have been provided concerning size 
distributions or other features of the topological structure of the vacuum. 

At the level of the Monte Carlo simulations, improved actions, which generally have 
better scaling properties than Wilson action, are expected to at least reduce the contribu¬ 
tion of dislocations to the susceptibility |^. Two proposals in this direction have been 
formulated recently. On one hand “hxed point perfect actions” have been constructed in 
terms of many loops and higher representations [^. These actions are dehned with the 
help of a blocking procedure whose choice should ensure the persistence of the good scaling 
qualities found at the fixed point even far from it, where the simulations are performed. 
On the other hand a proposal for a lattice action has been formulated using interpolating 
gauge fields, with the aim of ensuring on the lattice some of the properties of the continuum 
fields . Both of these approaches could be argued at the semi-classical level not to give 
a divergent topological susceptibility in the continuum limit. The first analyses of SU(2) 
topology using perfect actions have already been performed [^, |^. They indicate that 
even in this case there is an abundance of dislocations in the Monte Carlo conhgurations. 
To disentangle them from physically relevant structures a refined construction in many 
steps, using special reverse blocking procedures in combination with minimization steps 
has been developed to take advantage of the qualities of the hxed point actions |^, PB[ . 
By hxing the number of blockings (typically, one) one hxes the scale (in lattice units) at 
which one wants to extract the topological structure. 

We consider here the problem of measurement, although the lattice operators we use 
could also be used for the Monte Carlo generation of conhgurations, with great benehts. 
We propose to take advantage of the better scale-invariance of improved actions to design 
a cooling algorithm which by itself avoids most of the diseases of the usual cooling, and use 
it to analyze data from Monte Carlo simulations with any action. In fact, cooling should 
fulhll a number of requirements in order to qualify as an IR hltering procedure, that is, as 
a method to eliminate UV noise and permit studying topological excitations: 

- It should smooth out the short range huctuations, including “dislocations” 

- It should preserve the structure at all physical scales, including size and location of 
instantons under any amount of cooling 

- It should need no monitoring or any engineering which would only slow it down without 
ensuring stability and would make the algorithm conhguration dependent (introducing 
therefore a systematic error). 

The problem of finding a good cooling algorithm is of course related to that of finding an 
action possessing scale invariant classical solutions. From that point of view perfect actions 
provide a good basis for such an algorithm. However the use of many loops and higher 
representations makes cooling with the perfect actions computationally rather expensive. 
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Moreover the truncations performed to obtain an action easy to handle numerically may 
distort its good scaling properties (the truncated actions are no longer dehned rigorously by 
a hxed point equation). We look for a simpler construction of a good scale invariant action 
which does not involve as many loops as perfect actions and contains only linear dependence 
on the links, hence simplifying considerably the minimization algorithm. Since instantons 
are classical solutions tree level improvement is enough to guarantee scale invariance (up 
to corrections due to replacing integrals by discrete sums). We start therefore with order 
tree level improvement which is analytically dehned 1^, and account for higher order 
corrections by tuning a free parameter in the action to obtain instanton solutions which are 
practically scale invariant beyond some small-size threshold. This rather straightforward 
construction leads to an action with very good scaling properties (which therefore makes 
it interesting also as an improved action for Monte Carlo simulations). It also provides 
an improved topological charge operator. The improved cooling algorithm based on this 

2^) fulhlls most of the requirements stated 


(see also 28 


action and presented below 
above. The small size threshold is hxed in lattice units to ~ 2.3a (and therefore shrinks 
in approaching the continuum limit). It ensures that dislocations are eliminated during 
cooling while physically relevant instantons are preserved (including their size) if the lat¬ 
tice cut oh a is chosen small enough. Improved cooling is useful not only for producing 
good susceptibility data, but also for studying various conjectured ehects of instantons 
on the hadronic spectrum or on the chiral transition, by providing smooth conhgurations 
preserving the large scale structure of the original ones. The algorithm simply minimizes 
the action in the standard way and involves no further procedures (heating, blocking, etc). 

In the analysis concerning the determination of physical quantities like topological 
susceptibility, size distribution, etc we shall account for two potential problems, namely 


• threshold ehects (for coarse lattices 2.3a may already represent a physically relevant 
distance) and 


finite size ehects (small lattices may not accommodate large, physical instantons and 
moreover for certain sets of boundary conditions, self-dual Q = ±1 conhgurations do 
not exist in hnite volumes H). 


These questions will appear explicitly in our study of the scaling behavior and of the 
dependence on size and boundary conditions. To disentangle threshold from hnite size 
ehects we work at diherent values of the lattice spacing, a, with the same physical volume. 
We also use diherent lattice geometries (elongated and symmetrical) with periodic and 
twisted boundary conditions. This allows us to study threshold and hnite size ehects 
separately, and to assert that our lattices are such that the physical sizes of interest are 
ahected by neither of these ehects. 

Note that the only structures which in the end will survive cooling are those which are 
solutions of the lattice equations of motion. Instantons present in Monte Carlo conhgura¬ 
tions are only approximate solutions and during cooling will evolve toward true solutions. 
Some of the properties of the original instanton ensemble change therefore while cooling. 
These concern especially the features of the instanton - anti-instanton (TA) pair popula¬ 
tion (I-A pairs are not minima of the action, hence they annihilate). Since in the physical 
range of sizes practically no other changes in the topological structure occur during cooling 
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one can disentangle and stndy these effects to some extent by following the cooling his¬ 
tory. Thereby it is very helpfnl that we can nse the improved charge operator which allows 
to get meaningful charge density data already after a few cooling sweeps. The measured 
topological charge is an integer within 0(1%) already after about 5 sweeps. 

The paper is organized as follows: the improved cooling is introduced and analyzed 
in section In section ^ we present physical results for the topological susceptibility as 
well as for several properties of the instanton ensemble obtained by improved cooling on 
SU(2) configurations generated with the Wilson action. Both sections contain a number of 
necessary technical points, which however the reader interested only in the physical results 
may skip. Section ^ contains a summary of these results and our conclusions. 

While our analysis here refers explicitly to SU{2), most relations are valid for SU{N). 


2 The Method of Improved Cooling 

2.1 Improving the Action and the Observables 


As was pointed out in Ref. [27] the stability of instantons under cooling depends on the 


sign of the lattice corrections to the continuum action. On the ground of pure dimensional 
arguments it is easy to see that deviations from the continuum action for instantons of 
size p will go as S' — Stt^ = oq (a/p)^ + C>((a/p)"^), for p » a (from now on we will put a 
tilde over dimensionful quantities, otherwise everything will be given in lattice units, hence 
p = pa, etc). This allows to differentiate between “under-improved” (e.g., Wilson) actions, 
for which Oq < 0, and “over-improved” ones, Oq > 0; under the former instantons shrink. 


under the latter instantons beyond a certain size grow. This idea was used in to design 
a set of actions, S'(e), based on a combination of 1 x 1 and 2 x 2 plaquettes, with coefficients 

Clxl — — 7 ,— 5 C2x 2 — .o ’ 
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where e is a free parameter which controls the sign of the corrections. They are under¬ 
improving for e > 0 or over-improving for e < 0. Since these actions have been thoroughly 
studied in [^] we shall not discuss them here any longer. They have also been used in 


||3l| , |32| , to study topology on the lattice. 

Following this idea we addressed the problem of designing a lattice action for which the 
instanton action deviates as little as possible from the continuum one. Starting from the 
results in 0 one can construct a one-parameter set of actions with no O(a^) and 0{a‘^) 
corrections using hve fundamental, planar loops of size m x n: 

1 


x+nu 




mPri^ 






X 


x+mp 


S = 


( 2 ) 


2 = 1 


Here (m*, rii) = (1,1), (2, 2), (1, 2), (1, 3), (3, 3) for 7 = 1,..., 5 (other choices would imply 
larger loops) and: 


Cl = (19 - 55 C5)/9, C 2 = (1 - 64 c^)/Q 
C 3 = (-64 + 640 C5)/45, C 4 = 1/5 - 2 cg. 


(3) 
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Our choice for an improved action 


C5 = 1/20, 


(4) 


is obtained by tuning C 5 numerically, taking as criterion that the size of any instanton above 
a certain threshold po should remain unchanged during the cooling procedure (details will 
be given in section ^). We shall call this loop combination S{5Li), 5-loop-improved. We 


shall also discuss in section | 2 ]^ two further choices obtained by setting C 5 = 0 and C 5 = 1/10 
respectively, which we shall call “4Li” and “3Li”. 

In addition we use an improved topological charge operator based on the same set of 
five loops: 


Q = Ci 


i=l 


Q 


mi.rii- 


(5) 


with Ci as in eq. and where Qm,n is the naive topological charge operator defined through 
loops of size {m,n), i.e. 


Q 


m,n 


1 1 

327r^ 


(^Ff,i,{x]m,n)Fp^{x]m,n)'^ 


( 6 ) 


with Ffj_i,{x]m,n) given in terms of oriented clover averages of m x n plaquettes 


Fp^{x]m,n) 



(7) 


It can be easily shown that this definition of the topological charge is also free of 0(a'^) 
and 0{a‘^) corrections. We shall call Q{5Li) the combination for which C 5 is set as in eq. 
(1). The naive, non-improved, topological charge can be obtained from eq. by setting 
Cl = 1 and Cj^i = 0, we will call it in what follows Q{W). 

Our cooling algorithm exactly minimizes the local action at each step. It consists in an 
iterative replacement of each link by the normalized sum of staples connected to it through 
eq. (^. It involves no further calibration or engineering. The Q{5Li) charge operator 
produces reasonable density data already on still rough conhgurations and can be used to 
observe I-A pairs from the early stages of cooling. 


2.2 Instanton Description and Size Definitions 

Before describing how instantons evolve under our improved cooling it is useful to point 
out how the instanton size can be determined. This discussion provides also the basis for 
our analysis of the size distributions and other characteristics of the vacuum structure in 
the next section. 

One can design alternative definitions of the size based on the continuum’t Hooft ansatz 
1^, Q. This ansatz describes an isolated instanton in infinite volume. However, both the 
finite volume and the interaction among (anti-)instantons may lead to deviations from this 
ansatz (see, e.g. [^]). Moreover, the original “instantons” produced by Monte Carlo, 
which are only approximate solutions of the classical equations of motion, will typically 
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be deformed by quantum excitations showing up in a departure from spherical symmetry. 
To partially account for these effects in identifying these structures we use an ellipsoid 
generalization of the continuum ansatz for an instanton centered at {x°}: 


q{x) 


6N 


TT- 


n.=i Pm 


1 + 


4 

E 

M =1 


-I -4 


X,. 


X,. 


Pm 


( 8 ) 


and similarly for s{x) = s(a;)/87r^ (notice that we have introduced a normalized action S, 
and action density s(a;) dividing by Svr^). Here x^, are taken as real numbers and the 
ansatz eq. interpolates between the lattice points. For an instanton or anti-instanton 
the normalization iV is +1 or -1, respectively. The shape of the instanton is characterized 
concisely by the geometric average size 


p=n pF- 

11=1 

and the quantity: 



introduced as a measure of “excentricity”. 

We also use the partially integrated quantities (“prohles”): 


(9) 


( 10 ) 


Sm(^m) = E ^(^) 

(11) 



?m(2^m) = E ^(^)- 

(12) 


to be htted by the formula: 




3iV 

4Pm 


1 + 


C PV» PV »0 

Pm 


2 " 


- 5/2 




(13) 


Here periodicity effects are taken explicitly into account by adding in the £t satellites 
0^{x^ ± N^) {O = {q, s}), where is the lattice size. 

The size dehnitions we have used for our analysis are: 


(def.l) Ppeak, the so called “peak” radius, 


Ppeak (f?) 


6 


1/4 


7i‘^\q{x^)\ 


(14) 


with g(x°) the lattice topological charge density at the center of the instanton (one 
could use the action density s{x^) instead, we call the size extracted from it Ppeak(5')). 
To obtain the position of the center we interpolate between the lattice points using 
the ansatz eq. (H) (see (def.2) below). 
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(def.2) pc, extracted by solving eq. (P) on the central and nearest neighbor points for 
and then taking the geometric average eq. (|^). The corresponding ex- 
centricity (see eq. ([I0|) ) will be denoted e^- pc and Ppeak are related through eq. 

(S): 

Ppeak = (15) 


(def.3) pf, obtained with the help of a global fit of the configurations to eq. (H). For multi- 
instanton structures we use a linear superposition ansatz. Then we take again the 
geometric average eq. (j^) to obtain pf and we use eq. (p!0[) to obtain the excentricity 
ef. 


(def.4) Ppgak, based on partially integrated densities eqs. (^, |^, p!3|) : 

Ppeak(OM) = ( 16 ) 

with 0^(a;°) the value of s, q at the center of the instanton. Here we have interpolated 
between lattice points using a second order polynomial in x^. We call p^'' the size 
obtained from a full £t of the correspondingly integrated ansatz, eq. (0- 

(def.5) pd(.R)- For an isolated instanton, and assuming spherical symmetry as a hrst approx¬ 
imation, one can compute analytically the action density integrated over a hyper¬ 
sphere of radius R centered at the instanton position: 


(3i?2 + 

(p2 + R^f 


(17) 


This expression can be used to compute p given R and Sr. It will provide a set 
of definitions of the instanton size which we call pd(.R)- Of course for an isolated 
instanton well described by the continuum ansatz they should be R independent. 
The way we have proceeded to extract them from the data is as follows: locate the 
center of the instanton, compute the action contained in a sphere around it of radius 
R and use eq. ( p.2|) to obtain Pa{R)- The same can be done using the topological 
charge instead of the action. In particular we have used for our analysis R = Ppeak 
and R ppeak/ V^' 


All these dehnitions, most of which have already been used in other works, agree rea¬ 
sonably well with each other as long as the instanton does not depart too much from the 
continuum ansatz. Large instantons, which may be deformed by finite size effects, or very 
small instantons, distorted by lattice artifact effects, may give rise to differences. Of course 
most of the definitions above are based on a dilute instanton picture, since they do not 
take into account the overlap between instantons or instanton-anti-instanton interactions. 
Definition (def.3) based on a global fit to a linear superposition of continuum solutions eq. 
(1), can correct for the first of these effects as well as for periodicity effects, and will ideally 
account for all the topological objects present in the configuration. As a rule, after a few 
cooling sweeps the SLAdensities around the peaks of the action are well described by this 
ansatz. We will present below a detailed comparison of the different dehnitions for p, in 
connection with the Monte Carlo results. 






Of course, it is possible to use more “model-independent” definitions for the size, e.g. 
the 4-th root of the volume at whose boundary the density has dropped to half the value of 
the maximum (see, e.g., f^). However, the physical meaning and the compatibility of such 
dehnitions among themselves ultimately rely on the continuum ansatz. Since we found that 
the ellipsoid ansatz can take care of most cases we preferred to stay with dehnitions based 
directly on it in order to have less arbitrariness in the interpretation of the results and to 
fully automatize the analysis (with exception of (def.3) all size determinations above can 
in fact be done on line). 


2.3 Tuning of the Action for Optimal Scale Invariance 

To check the effects of tuning the parameter C 5 in eq. (^, we cooled with different values 
of C 5 a set of instantons of various sizes (generated numerically by cooling with ^(e) with 
various e in order to vary their sizes), obtained on 12^ lattices with gauge group SU{2). 
Since we want to disentangle small distance from hnite lattice size effects (see next section), 
we use here twisted boundary conditions in the time direction with k = ( 1 , 1 , 1 ) twist, to 
ensure that instanton solutions exist on hnite lattices PU, 


In Fig. ^ we compare the 
cooling behavior for diherent choices of C 5 of: (a) the action and (b) the size Ppeak (def.4) for 
the same instanton, as a function of the cooling step. The two improved actions denoted 
3Li and 4Lz correspond to the settings C 5 = 1/10 and C 5 = 0, respectively, in eqs. (|^, 
^). Although at tree level both of them have no 0(a‘^) and Cl(a^) corrections they tend 
to act over- or under-improving, respectively, while the choice bLi preserves the size of 
all instantons above the threshold for an indehnite number of sweeps. This illustrates the 
signihcance of the tuning eq. (^. The 3Li loop combination, however, is simpler and since 
it also guarantees that physical instantons will not decay under any number of cooling 
sweeps, it may be useful for more extensive analyses (large lattices, etc). 

Fig. 1^ shows the behavior of instantons of various initial sizes under SLi-cooling. They 
remain practically unchanged over any practicable number of cooling sweeps provided that 


p> po^ 2.3, 


( 18 ) 


(in units of a). The size of the original instanton has been varied by applying S{e) cooling, 
therefore the original instantons do not correspond to minima of S{5Li). Submitting such 
conhgurations to bLi cooling hrst re-adapts them to the new equations of motion, implying 
the small changes which can be observed in the hrst ~ 50 sweeps of Fig. ^ 

In Fig. |] we present a comparison between the size dependence of the 5LTaction and 
Wilson action. Instantons of different sizes are extracted from various stages of Wilson 
cooling of 3 different starting conhgurations. The triangles are data from Monte Carlo 
instantons (/? = 2.5, 12^ lattice) cooled with 300 sweeps of bLi cooling. For the improved 
action the dependence on the size is close to a step-function, giving to better than 
0 . 1 % for p > po and a steep descent (slightly conhguration dependent) in the interval 
(0.8 — 1.0)po |@- As Fig. ^ shows clearly, there is a small energy barrier which prevents 
the decay of any instanton of size greater than po- 

There are simpler actions with no O(a^) corrections, like the e action described in 


section | 2 . 1 | with e = 0 , which have been shown to improve results on the susceptibility 
(see for a QCD study using such an action). Our analysis indicates, however, that the 
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(9(a^) terms may be relatively important in the region of physical sizes and they may lead 


to shifts in the instanton size with cooling, as already observed by Brower et al [36|. These 


terms are also important in the threshold region, which tends to loose its step-fnnction like 
behavior. The cooling dependence of the sizes will also make a scaling test more involved. 


2.4 Finite Size Effects 


Lattice instantons may differ from the continunm ansatz due also to the hnite lattice 
volume, therefore the question of stability of instantons cannot be completely decoupled 
from hnite size effects. The inhuence of hnite size ehects on the quantities we want to 
study comes mostly in two ways: 


They ahect the stability of instantons under cooling. Since on a periodic hyper-torus 
self-dual Q = ±1 solutions do not exist H. single instanton conhgurations will 
always collapse after a while during cooling. Of course if the lattice size is much 
larger than the instanton size there is a wide window of meta-stability for Q = ±1 
which for all practical purposes would amount to stable instantons. This window 
can be enlarged by choosing at least one of the lattice sizes, say Nt, very large. The 
problem can be completely avoided by imposing twisted boundary conditions in one 
of the lattice directions. In this case our cooling algorithm guarantees stability of the 
cooled structures for any value of the topological charge and an indehnite number of 
cooling sweeps. 


(b) There is a maximal instanton size Pmax, determined by the length of the lattice and 
by the boundary conditions, hence the large-size part of the size distribution will 
be cut off and distorted. Due to point (a) above, pmax will also depend on the 
total topological charge of the conhguration. We have observed that for self-dual 
Q = ±1 conhgurations and for symmetric lattices with periodic boundary conditions 
(p.b.c.) p^ax ^ Ns/5 1^, for elongated lattices with p.b.c. p^ax ^ Ng/A and for 
twisted boundary conditions (t.b.c.) p^ax ~ 2Ns/5. Moreover large instantons which 
“feel” the boundaries will typically be deformed by them. As a consequence for large 
instantons different size dehnitions can give quite different results. To obtain a size 
distribution free from this effect the typical physical instanton size has to be less than 
the maximal instanton size, pmax, allowed by the lattice. Note however that for higher 
topological charge pmax can be larger and may not show such a strong dependence 
on boundary conditions, hence the previous estimate is quite conservative. 


This provides a general window to retrieve the physical size distribution, which goes 
from the cooling threshold po = 2.3a to at most the maximal instanton size. The threshold 
is a feature of the algorithm and is essential in discarding dislocations; it is fixed by choosing 
(3. As already said pmax is fixed by the lattice size, boundary conditions and topological 
charge of the configuration. 

Point (a) is illustrated for 12“^ and 12^ x 36 lattices on Fig. ^ where we also show 
for comparison the behavior of a Q = 2 configuration. For Ng = Nt = 12 p.b.c. the 
maximal instanton size (attained by over-improved cooling) is too close to our stability 
threshold po = 2.3, in consequence any isolated instanton decays quite fast. We expect 
the region of meta-stability to be enlarged when pmax gets away from the threshold (see 
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for instance the 12^ x 36 p.b.c. lattice). We stress that the instability only affects self- 
dnal Q = ±1 conhgnrations bnt not mnlti-instanton conhgnrations and can be completely 
avoided by nsing twisted bonndary conditions in at least one lattice direction. The effect 
on the snsceptibility and size distribntion measnred on p.b.c. lattices will be less strong for 
large physical volnmes where the contribntion to the snsceptibility of conhgnrations with 
single instantons is smaller. 

As an illnstration of the nncertainties in assigning physical parameters to very large 
instantons we present in Table |2.4| data coming from the analysis of the maximal size 
instanton on a 12^ lattice with twisted bonndary condition in time. Althongh this (com¬ 
pletely stable) conhgnration has total charge 1 and action Stt^ the nnmber of peaks in 
the energy density is 6 and in conseqnence a naive assignment of instantons to peaks will 
misleadingly interpret it as 6 different instantons, with qnite large sizes and excentricity. 
On the other hand an analysis of the time prohle reveals only one peak with mnch smaller 
size. 


{x,y,z,t) 

Pc 

Pp 

Cc 

N 

(3,4, 2, 2) 

7.90 

7.22 

0.66 

-1.43 

(4, 6, 2, 2) 

7.61 

7.09 

0.58 

-1.33 

(1,1, 3, 2) 

7.38 

6.96 

0.58 

-1.26 

(6, 7, 3, 2) 

7.31 

7.00 

0.45 

-1.19 

(10,4,5,2) 

10.06 

7.52 

2.21 

-3.21 

(10,10,5,2) 

7.93 

7.23 

0.66 

-1.44 


Table |2.4| : Assignment of parameters made by the htting programs to the maximal size 
instanton {S = 1.0000, Q = —1.0001) on a 12"^ t.b.c. lattice. {x,y,z,t) represents the 
location of the 6 peaks in the action and charge densities. N is the normalization appearing 
in the continnnm ansatz, eq. (§), as given by the “center and nearest neighbors” £t (def.2) 
in section |2.2| . Notice that all the peaks have the same time coordinate, hence only one 
peak will be detected in the integrated densities St{t) and qt{t) eq. (|^ giving a valne 
Ppeak = 4.66 (the time coordinate is selected by the twist). 


3 SU(2) Topology by Improved Cooling 

3.1 Simulation with Wilson Action 

The hrst criterion for the physical relevance of all resnlts is how they scale with the re¬ 
moval of the cnt-off. We work in a region of (3 where one expects to see scaling behavior, 
therefore all qnantities of physical signihcance (snsceptibility, charge distribntion, size dis¬ 
tribntion, etc) shonld scale according to their dimension when approaching the continnnm. 
As mentioned above, one may expect short-range topological excitations which are nn- 
physical in the sense that they do not correspond to attribntes of the continnnm physics. 
These are therefore artifacts of the discretization and will show np as non-scaling, cnt- 
off dependent strnctnres in the resnlts. Onr algorithm has an inherent cnt-off for snch 
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structure, the threshold po — 2.3. Since the threshold is hxed in lattice units, it shrinks 
physically and it should not affect physical sizes if a is small enough: the physical threshold 
Po ~ 2.3a —>■ 0 for a —> 0. This we shall test below when studying the scaling of various 
quantities. On the other hand, since Monte Carlo simulations with the Wilson action copi¬ 
ously produce short range fluctuations one can generally ask whether this shrinking cut-off 
at Po can indeed ensure the absence of unphysical fluctuations in cooled configurations 
when approaching the continuum. 

Following an argument of Pugh and Teper H we write the contribution of small 
instantons of size po to the partition function as (A^ = 2 here) 

[po(a)a]-'e-^-("°)/^^ a{(5) ^ A’^e'^, A = (19) 

The contribution of instantons of a size p such that Sy/{p) < l/(27r^A^/3i) = 24/llA^^ 
diverges for f3 = Ajg^ —^ oo. We can calculate Sw{p) by cooling large instantons with 
the Wilson action and measuring the action and the size as the instantons shrink during 
cooling (see Fig. 3). This analysis leads to an UV threshold 1 < po < 2. It is difficult to 
get a precise estimate since the dependence of Sw on p in this region is sensitive to details 
of the configuration. We can conclude, however, that our threshold po of about 2.3 is fully 
satisfactory when dealing with configurations produced with the Wilson action: it ensures 
removal of the divergent contribution of lattice artifacts, without cutting too much from 
the physical sizes even at rather low (3. 

Likewise, the correction due to the finite lattice size should be estimated by changing 
the volume and/or boundary conditions and checking the sensitivity of the results to these 
changes. Quite generally, they will affect structures at the long-range edge of the size 
distribution. To estimate the importance of the finite lattice size effects we use both 
different physical volumes and boundary conditions. 

We base our analysis on the following SU{2) Monte Carlo simulations (heat-bath, 
Wilson action): 



Ns 

Nt 

/? 

a{j3) (fm) 

L (fm) 

b. cond. 

#conf. 

#sw. betw. conf. 

(la) 

12 

12 

2.4 

0.12 

1.44 

t.b.c. 

237 

100 

(lb) 

12 

12 

2.4 

0.12 

1.44 

p.b.c. 

202 

100 

(Ic) 

12 

36 

2.4 

0.12 

1.44 

p.b.c. 

198 

100 

(2) 

12 

12 

2.5 

0.085 

1.02 

t.b.c. 

160 

250 

(3) 

24 

24 

2.6 

0.06 

1.44 

t.b.c. 

84 

200 


Table Ensembles of gauge field configurations analyzed; t.b.c lattices have k = (1,1,1) 


twist in the time direction and periodic boundary conditions in space; p.b.c. lattices have 
periodic boundary conditions in all directions. 


where we have estimated a(/9) using the 1-loop phenomenological formula (derived from 
data in [^ ) 

a(^()=400exp|-^}fm. (20) 
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Here all data are taken after 20000 sweeps for thermalization. Part of these data have been 
presented at the Lattice’95 and Lattice’96 conferences. 

As a rule, during cooling we performed measurements after 5, 20, 50, 150 and 300 
cooling sweeps for lattices (1) and (3), and every 20 sweeps for lattice (2) (up to 300). 

We typically calculate the bLi charge and electric and magnetic action 4-dimensional 
densities and thereby the total topological charge and action of the conhgurations. 


3.2 Topological Susceptibility and Charge Distribution 

Since the topological charge Q{5Li) stabilizes very fast to an integer within less than 
1 %, we expect deviations from the physical values in the susceptibility to be only due to 
threshold or to finite size effects. Our results are presented in Fig. ^ and Table |3.2| a. The 
susceptibility (from Q{5Li)) settles early in the cooling. It scales very well with the cut¬ 
off (compare the (1) and (3) data) and shows little dependence on the volume (compare 
lattices (la-b) and (3) with (Ic) and with (2)). Data coming from the smaller physical 
volume seem to be higher by about 5%, although they are compatible within errors with 
the other two determinations. The effect can be due to our particular estimation of the 


physical scale eq. (|20|) or it can be a volume effect (that the lattice size for these data 
is too small is evident in the size distribution - see next section). The small decrease 
in the susceptibility during cooling showing up in the data from lattices (1) (/? = 2.4) 
appears to be mostly a threshold effect. Our cutoff on small sizes (po = 2.3a) will induce 
a systematic underestimate in the measurement of the susceptibility since some physical 
instantons might be eliminated, but this error goes down with increasing and can be 
easily estimated using the size distribution. As we will see in section |3.3| we obtain size 
distributions which scale nicely. Based on the £t to these distributions this underestimate 
should be about 3% for a = 0.12fm and less than 1% for a = 0.06fm. This is compatible 
with the difference between short and long cooling values (see Table |3.2| a.) and does not 
change signihcantly the overall picture. Assuming the Witten - Veneziano formula to hold, 
our result = 200(15)MeV agrees excellently with the phenomenological expectation. 

Recent results for SU(2) are about 20% higher (for x itself they are about 

a factor of two larger). We believe the discrepancy might be caused by the presence of 
some dislocations in the data of [25|. It would be very interesting to have data for the 


size distribution extracted using inverse blocking. A scaling check on such distribution will 
provide an unambiguous way of deciding whether dislocations are still present, since their 
contribution to the size distribution does not scale. 

Further confidence in the correctness of our results for the topological susceptibility has 
been provided by R. Narayanan and P. Vranas [|^. They have used the overlap formalism 
m to measure the topological susceptibility on a 12^ lattice at /3 = 2.4 with periodic 


boundary conditions (our (lb) lattice). Their results for the topological susceptibility as 
well as for the charge distribution are in perfect agreement, within errors, with ours. 

Note that the early stabilization of the susceptibility with cooling is not only due to 
the efficiency of the improved algorithm in removing dislocations, but also to the fact 
that we use the improved charge operator Q{5Li). To illustrate this we present in Ta¬ 
ble |3.2|b the values of the topological susceptibility, for lattices (la) and (2) in Table O 


extracted from the measurement of Q{W) over the same set of bLi cooled conhgurations. 
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(la) 

(lb) 

(Ic) 

(3) 


(2) 

sw 

5Li 

5 Li 

bLi 

5Li 

sw 

bLi 

0 

161(5) 




0 

219(8) 

5 

199(3) 

200(5) 

201(5) 

205(16) 



20 

200(3) 

200 (6) 

199(5) 

199(15) 

20 

212(7) 

50 

198(3) 

198(6) 

195(5) 

197(15) 

60 

212(7) 

150 

196(4) 

195(6) 

192(6) 

197(15) 

140 

212(7) 

300 

195(4) 

193(6) 

189(6) 

194(14) 

300 

211(7) 


Table 3.2| a: (MeV) computed from the improved charge operator Q{5Li), vs number 

of cooling sweeps. Errors have been estimated using the jackknife method. 


The susceptibility stabilizes much earlier with cooling if the improved operator is used. 
Indeed even when no cooling at all is performed, the deviation from the asymptotic value 
is at most 20% for Q{5Li) while it is twice as much for the naive charge operator (after 
5 sweeps Q{5Li) is already asymptotic while Q(W) is still 10% off). The results obtained 
with Q{5Li) after only a few cooling sweeps deviate by no more than 0.5% from what is 
obtained by approximating the charge by the closest integer. Even without cooling the 
deviation is only about 1 — 2%. This rounding procedure also works for Q(W), but only 
after a long cooling: approximating the charge by the closest integer after 300 sweeps we 
obtain = 194(4), 212(7) (MeV) for (la) and (2) respectively, in perfect agreement with 
the results from Q{5Li). The importance of improving the topological charge operator for 
non-smooth configurations has been recently stressed in . 



(la) 


(2) 

sw 

Wilson 

sw 

Wilson 

0 

123(3) 

0 

167(6) 

5 

179(3) 



20 

186(3) 

20 

201 (6) 

50 

188(3) 

60 

204(7) 

150 

188(4) 

140 

204(6) 

300 

187(4) 

300 

204(6) 


Table p.2| b: (MeV) computed from the naive charge operator Q{W) vs number of 

cooling sweeps. 

The topological charge distribution is presented in Fig. for various boundary condi¬ 
tions and numbers of cooling sweeps. The differences induced by the former are related 
mostly to the instability of self-dual charge ±1 conhgurations for periodic boundary con- 
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ditions (the p.b.c. data show typically more conhgurations at Q = 0 and fewer at Q = ±1 
than the t.b.c. data, bnt the statistics is not snfficient for a signihcant signal). The dis¬ 
appearance of some narrow peaks (beyond pair annihilation) when going from 20 to 300 
cooling sweeps does not lead to significant changes in the distribntion. The shape of the 
charge distribntion agrees well with a ganssian of width given by the snsceptibility. These 
data are also nicely corroborated by those from |T| . 


3.3 Instanton Shape and Size Distribution 

For a systematic analysis of the local structnre we developed a nnmber of programs to 
automatically recognize and evaluate the local features. We think it useful to give some 
technical details about the procedure, to permit an overview of the various constraints and 
uncertainties. 

Essentially these programs proceed by Ending first the peaks in the action and in the 
topological charge density. The identification of the peaks proceeds along the following 
three steps: 

• Look for peaks in the electric and magnetic parts of the action density, E‘^{x) and 
B‘^{x). Assign a peak to every lattice point whose density is maximal with respect 
to its nearest neighbors. We also assign a normalization ±1 to each peak depending 
on whether its topological charge is positive or negative. 

• Apply a cut that removes from the set of peaks those which are so fiat that Ppeak(5') 
is larger than Ng. Typically small fluctuations, not associated to instantons, or 
remnants of I-A annihilation will give rise to very low peaks in the energy density, 
which misleadingly translate via the Ppeak formula eq. into very large sizes. With 
this cut we suppress part of these effects. 

• Apply a cut from self-duality by comparing the locations of peaks in E‘^{x) and B‘^{x) 
(with equal normalization). Only peaks for which the distance between locations is 
less than or equal to 2a are kept. 

After performing these cuts we have observed that there is an abundance of close peaks 
separated from each other no more than 1 lattice spacing along each direction. This be¬ 
havior has been also pointed out in . We have adopted the attitude of identifying such 
peaks as being only one. We believe that the effect is mostly due to locating peaks by 
looking only at nearest neighbors and not at points separated along the diagonal (slightly 
excentric instantons located diagonally between the lattice points may lead to such a sig¬ 
nal). The effect of this last cut reduces the number of peaks by about 8% independently of 
the number of cooling sweeps. It is important to point out that the cut is applied in lattice 
units so if there is indeed a physical effect associated to close peaks it should be possible 
to detect it anyway by reducing the lattice spacing. In section p.4.2| we will discuss further 
this effect by comparing results from lattices with a = 0.12 fm and a = 0.06 fm. 

The typical level of accuracy achieved by these programs turned out to be above 90%, 
the failures being mostly related to structures affected by threshold or by finite size effects 
(e.g., wrong normalization, misinterpretation of a wide, deformed instanton as a multiple 
peak structure, etc). While global quantities like total charge or topological susceptibility 
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are not affected by such uncertainties, the latter will be reflected in systematic errors 
concerning, e.g., the small and large size edges of the size distributions. The errors we 
quote in this paper are statistical only. But we now proceed to estimate the systematic 
errors. 


3.3.1 Determination of Instanton Shape and Size 


Both hnite size effects and instanton - (anti-)instanton interaction may deform the spher¬ 
ically symmetric, continuum ansatz and we did not want to bias a priori the information 
about instanton shape by assuming spherical symmetry. We decided therefore to evaluate 
simultaneously several of the size dehnitions in section | 2 . 2 | . A comparison between them 


and more specihcally a measurement of the excentricity (see eq. ([T0|)) provides a good 
estimate of the departure from the continuum ansatz. Our analysis proceeds along two 
directions: 

(a) Extracting “on line” sizes corresponding to dehnitions 1, 2, 5 in section (i.e. Ppeak, 
pc and pd{R) for R = Ppeak and R = Ppenk/V^)- We also extract the excentricity €c 
using (def. 2 ) in section ^172. 


(b) Performing by MINUIT global hts of the bLi charge and action densities using eq. (§). 
In most cases this ansatz hts rather well the SLAdensities (charge and action) of 
isolated instantons already after a few cooling sweeps. The simple superposition 
extension of this ansatz can also handle overlapping structures in multi-peak hts 
over clusters. 

We tried to assess the uncertainty of the analysis self-consistently by comparing all the size 
dehnitions. 

The global hts (b) start with initial information about the location and size of the 
peaks provided by (a). Clusters are hrst identihed, by considering as neighbors peaks at 
a distance less than their average radius. For this we use the size estimate from (a), or 
a fraction thereof if we want to limit the maximal cluster size (MINUIT would not cope 
with more than 5 peaks at a time). Once identihed, the clusters are htted separately with 
a superposition ansatz. The MINUIT ht uses 8 or 9 parameters per peak (in the latter 
case the normalization N in eq. (H) is taken as a real number and permitted to vary). 
The htting procedure tends to fail in the cases where also the peak assignment is not well 
dehned (see the discussion in sect. |^) . For the description of the topological structure 
uncovered in this way we shall use the geometric average size eq. (j^) and the “excentricity” 
eq. ( 0 . Since the global hts are rather expensive we performed them only on subsets of 
10 conhgurations for the (la) and (lb) lattices and 5 conhgurations on (Ic) (for the 24"^ 
lattices we could not perform global hts on our workstations), therefore the data for pf 
and 6 / ((def.3), section ^(^) have large statistical errors. 

A comparison between the ht analysis, Ppeak, (def.l) and pc (def.2) is given in Table 
|3.3.1| a. The diherences in the results for size and the excentricities are presented as the 
root mean squares V/p|rms and e/|rms, with 5p/p = 2(ppeak -p)/(Ppeak + p)- Here p = pc or 
p = Pf as obtained from the following hts: “f4” (free norm cluster ht with up to 4 peaks in a 
cluster), “fn” (cluster ht with up to 4 peaks in a cluster with the norm N = ±1 as assigned 
by (a)) and “fs” (single peak ht with the norm N taken as free parameter). The hgures in 
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brackets correspond to discarding peaks for which one or more of the sizes are larger than 
the lattice size - these cases typically prove to be I-A pairs in the course of annihilation or 
simply failures of the £t program in the presence of very distorted structures. They affect 
5 — 10% of the peaks. 

As can be seen, the fitted radius averaged over directions agrees quite well with Ppeak 
(def.l). From the data in Table |3.3.1| a we estimate a total error of less than 10% for long 
cooling: we reproduce this (informal) estimation a{p) in Table 3.3.1 b. Both (j{p) and 
e/|rms decrease with the number of cooling sweeps, which is understandable as an effect of 
smoothing and of the disappearance of I-A pairs. 

In Fig. 1^ we compare the size distributions obtained from the Ppeak and Pc data. We 
also show here the “excentricity” distribution Cc. As is evident in the cross-correlation 
between size and excentricity on Fig. |^, larger instantons depart more from the spherical 
shape. Also, the average excentricity decreases with cooling. The excentricity can have at 
least four origins: 


• sensitivity of the £t procedure to details of the conhguration not considered in the 
ansatz, e.g., the presence of nearby instantons; 

• interaction among (anti-)instantons, which is not taken care of by the ansatz; 

• physical fluctuations (e.g., higher momentum states) in the ensemble of instantons 
produced in the Monte Carlo simulation at physical sizes; 

• distortion of the true lattice instantons at large sizes due to hnite volume effects. 


Both the hrst and second effect are stronger if the instantons have a large overlap, therefore 
they should decrease with cooling as a result of the dilution of the ensemble due to pair 
annihilation. For the same reason they also could affect large instantons more than small 
ones. The smoothing out of physical fluctuations during cooling will tend to decrease the 
average excentricity for instantons in the range of sizes not affected by hnite volume effects, 
since there the exact solution is spherically symmetric; this does not necessarily hold for 
large instantons, which may be very deformed asymptotically. Fig. shows that excentric 
instantons are rare. We have explicitly verihed that the 4 axes p^ of our ellipsoidal hts are 
strongly correlated. Departure from spherical symmetry is moderate, and well accounted 
for by our ellipsoidal ansatz, eq.(|). 

In Fig. ^ we compare sizes obtained from dehnitions 1 and 5 (see section |^) : Ppeak(*S'), 
Ppeak(Q) or pd{R)- All should agree if the continuum ansatz is correct. As we can see, the 
differences are typically well within the uncertainty estimated on the basis of the hts, and 
mostly ahect large instantons of size above Ng/2. The ehects related to large instantons 
are clearly rehected in Fig. ^ where we compare size distributions after 300 cooling sweeps 
extracted from Ppeak(Q) (def.l) and p\ (def.4) for lattice (2) (12^ at /3 = 2.5, physical size 
approximately Ifm). The Pf distribution has a strong large size cutoh around Ns/2 while 
the Ppeak (Q) one is quite biased towards large sizes. In addition the number of peaks in the 
energy density is larger by as much as 40% than the one extracted from spatially integrated 
densities (it decreases from 98 to 61 after integration for 70 configurations analyzed). 

A further measurement of the good agreement between diherent size dehnitions is pro¬ 
vided by the value of the normalization entering in the density formula eq. (8) (see also eq. 
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(15)). In Fig. 1^ we present the histogram for the distribution of normalizations obtained 
from the “center and nearest neighbors” analysis (def.2). It is nicely centered about ±1 
with a dispersion that decreases with the number of cooling sweeps due to effects analogous 
to the ones quoted above in relation with the excentricity. 


sw 

£t 

(la) t.b.c 
#peaks ^Irms 

f/ rms 

(lb) 

Speaks 

p.b.c. 

5p 1 

p |rms 

rms 

(Ic) 

Speaks 

p.b.c. elongated 

^1 eA 

p |rms |rms 


c 

98 

0.125 

0.427 

103 

0.117 

0.333 

172 

0.079 

0.552 



(96) 

(0.115) 

(0.318) 

(101) 

(0.109) 

(0.290) 

(171) 

(0.076) 

(0.552) 


fn 

98 

0.125 

1.427 

103 

0.203 

1.563 

172 

0.099 

0.622 

20 


(96) 

(0.123) 

(0.386) 

(94) 

(0.205) 

(0.375) 

(168) 

(0.093) 

(0.414) 



98 

0.176 

0.536 

103 

0.191 

0.962 

172 

0.184 

0.808 



(95) 

(0.170) 

(0.326) 

(101) 

(0.158) 

(0.308) 

(165) 

(0.176) 

(0.404) 


fs 

98 

0.140 

0.527 

103 

0.157 

0.352 

172 

0.144 

0.384 



(97) 

(0.136) 

(0.293) 

(101) 

(0.158) 

(0.308) 

(171) 

(0.143) 

(0.375) 


c 

17 

0.033 

0.211 

15 

0.246 

0.248 

28 

0.065 

0.221 



(17) 



(13) 

(0.084) 

(0.115) 

(28) 




fn 

17 

0.095 

0.295 

15 

0.222 

0.253 

28 

0.100 

0.409 

300 


(17) 



(13) 

(0.083) 

(0.153) 

(27) 

(0.060) 

(0.211) 



17 

0.101 

0.275 

15 

0.179 

0.274 

28 

0.064 

0.246 



(17) 



(13) 

(0.021) 

(0.127) 

(28) 




fs 

17 

0.101 

0.272 

15 

0.238 

0.391 

28 

0.092 

0.274 



(17) 



(13) 

(0.021) 

(0.121) 

(28) 




Table |3.3.1| a: Comparison of £t sizes and excentricities 


sw 

(la) 

^(P) 

(lb) 

(Ic) 

(^{p) 

20 

300 

0.14 (0.12) 
0.06 (0.06) 

0.16 (0.15) 
0.20 (0.05) 

0.14 (0.14) 
0.08 (0.06) 


Table 3.3.1b: Informal estimation of systematic errors 
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A general conclusion from the analysis described above is the good agreement, up to 
at most 10% discrepancy, of all the size definitions. Discrepancies appear due mostly to 
hnite size effects and I-A pairs in the late stages of annihilation, which typically reflect in 
a departure from spherical symmetry and a wrong assignment of peaks. Based on these 
observations, we decided to quote in the following results for Ppeak (def.l, eq. (|T^) only. 
(In 1^ the instanton sizes for lattice (2) were extracted from the t-profiles - see eq. (p!6D). 
Systematic errors in the determination of the size are estimated to be at most 10%. 

For the lattice (3), we used lattice coordinates in eq. (0) to define the instanton center 


and no continuous interpolation of 


in Xp was attempted. 


This may lead to an 


overestimation of the instanton size by up to 1/p^ which for the smallest stable instantons 
allowed by our algorithm, p = 2.3, amounts to at most about 15%, while for the typical 
size instanton it is below 2%. 


3.3.2 Instanton Size Distribution 

The results for the size distributions are presented in Fig. We only compare the 
size distributions for lattices (1) and (3) since for lattice (2) the large-size edge of the 
distribution is strongly affected by finite size effects (see discussion in section |3.3.1| and 
Fig. P{p) counts the number of instantons at distance p normalized in such a way that 
the integral over p is equal to 1. 

The primary criterion for the physical significance of these distributions is how they 
scale in approaching the continuum limit. As we see from Fig. 0 the threshold (indicated 
by the vertical dashed line) is gradually removed with decreasing a and the distributions 
scale nicely above it. The data at /? = 2.6 show that most of the physically relevant 
structure is in a well dehned region above the threshold, where the latter has no influence. 
Comparison of the (1) and (3) data shows that at /? = 2.4 the bulk of the size distribution 
already begins to get out of the threshold region. The slight decay of the susceptibility 
with cooling observed in Table |3.2| a and Fig. ^ for (3 = 2.4 can be completely understood 
as due to the disappearance of a few, small but still physically relevant instantons close 
to the threshold (compare in Fig. 0 the distribution of instantons close to the threshold 
between 20 and 300 cooling sweeps). This good scaling of the distributions when changing 
the lattice cut off by a factor 2 gives evidence for their physical relevance. The stability 
under cooling of the distributions is also a check of the physical relevance of the results 
and the success of our algorithm in preserving the structure of the original conhgurations. 


In Table |3.3.2| we give the data for the non-normalized size distributions, the value pM 
of the instanton size at the peak of the distribution and the “width” w of the distribution 
extracted from a £t 

F(p)oc (21) 

which for small p reproduces the semi-classical result 


p(p) oc y-'® 


( 22 ) 


and for large sizes introduces a damping which cuts the infrared divergence. This ansatz 
can be derived by using the semi-classical distribution with a modihed running coupling 
constant which stops running beyond some value of p (see P). Other plausible ansatze, 
like a Gaussian distribution, or p'^P times a Gaussian, give significantly poorer fits. 
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The curve in Fig. 0 represents a fit with eq. (^if) to the average of all (1) data after 
20 cooling sweeps and the (3) data after 50 cooling sweeps (leaving out the data below 
threshold for (1)) with a per degree of freedom of 0.9. This £t gives a typical instanton 
size pM = 0.43(5) fm, a full “width” w = 0.47(9) fm and a value p = 3(1) for the exponent. 

If the action used in the Monte Carlo simulation introduces spurious structure above the 
threshold, this structure will of course be seen by our algorithm, and should be identified 
as such by a scaling analysis. Since the Wilson action S'(hF) for instantons of size po is 
lower than 1 (see Fig. ^ we may expect it to produce in the Monte Carlo simulation an 
unwanted excess of instantons there. This may lead to an “artificial” increase in the size 
distribution in the region just above Pq. This effect should disappear when the physical 
sizes involved move away from the threshold. From Fig. 0 we see that the relevant part 
of the size distribution scales. This is indication that the effect under discussion is small, 
and that the size distribution we measure may be considered physical. 


p (fm) 

0.18 

0.30 

0.42 

0.54 

0.66 

0.78 

0.90 

1.02 

1.14 

1.26 

1.38 


Pm 

xVdof 

(la) 20sw 

0.55 

2.50 

3.24 

2.80 

1.14 

0.34 

0.01 

0.00 

0.00 

0.00 

0.00 


0.43 

1.0 

300sw 

0.00 

0.36 

0.48 

0.40 

0.37 

0.21 

0.17 

0.04 

0.04 

0.03 

0.01 


0.44 

0.6 

(lb) 20sw 

0.60 

2.46 

3.37 

2.72 

1.18 

0.23 

0.01 

0.00 

0.00 

0.00 

0.00 


0.43 

0.3 

300sw 

0.03 

0.47 

0.68 

0.31 

0.16 

0.07 

0.01 

0.01 

0.01 

0.00 

0.00 


0.38 

1.1 

(Ic) 20sw 

1.82 

7.85 

9.48 

7.87 

3.70 

0.80 

0.05 

0.00 

0.00 

0.00 

0.00 


0.43 

1.6 

300sw 

0.07 

0.73 

1.18 

1.16 

0.70 

0.29 

0.17 

0.07 

0.06 

0.04 

0.04 


0.46 

3. 

(3) 50sw 

2.36 

3.26 

4.32 

3.65 

1.58 

0.44 

0.08 

0.02 

0.01 

0.00 

0.00 


0.39 

2.5 

300sw 

0.57 

0.98 

0.82 

0.76 

0.40 

0.25 

0.12 

0.02 

0.01 

0.00 

0.00 


0.37 

0.5 


Table |T3^ : Non-normalized size distributions and typical instanton size (we only divide by 
the number of configurations; notice that we have changed the binning for the 24"^ lattice 
with respect to Fig. O). For the bins p = 0.06fm and below there are no contributions. 


The typical instanton size we obtain can be compared with previous results obtained in 
Refs. and more recently and Polikarpov and Veselov cool using the 

Wilson action for a small number of sweeps and quote an average size of about 0.38 fm. 
Michael and Spencer use a cooling algorithm based on the Wilson action and adjusted 
to provide a meta-stability window for the instantons. They have results on 16^ and 24^ 
lattices at /? = 2.4 and [3 = 2.5 respectively. The peak of their size distribution depends on 
the lattice size, obtaining pM = 0.38 fm and 0.26 fm for 16^ and 24^ respectively. We believe 
that their method is subject to systematic errors induced by the unavoidable shrinking of 
instantons under Wilson cooling. They have carefully adjusted their algorithm to obtain 
the same “degree” of cooling on both lattices but still the difference is large. Our method 
provides a safe way of ensuring no distortion of physical structures without need for any 
further tuning or calibration of the algorithm. Chu et al use a convolution with the 
continuum ansatz of the density correlations obtained by cooling with Wilson action for 
the SU(3) theory. They quote an instanton size of about 0.36 fm after 20 cooling sweeps. 
Brower et al use an 0{a‘^) improved action in connection with a relaxation variant of 
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cooling for SU{3) theory without and with dynamical fermions. They apply an instanton 
identihcation procedure also based on the continuum ansatz similar to the method used 
~9| and further developed here. Since the maximum in the size distributions shown in 
611 drifts from about 0.4 fm after 20 “cooling” sweeps to about 0.5 fm after 50 sweeps 


m 


(presumably as a result of the 0{a^) terms in the action) the comparison is however difficult. 

Finite size effects can be estimated by comparing (la,b,c). The results agree perfectly 
for short cooling. After 300 cooling sweeps there are several sources of systematic errors 
which can give rise to slight differences. For lattices (la,b,c), a = 0.12fm, and there are 
some physically relevant instantons close to the threshold which disappear with cooling. 
This fact tends to induce a small overall shift of the (1) size distributions towards larger 
sizes. This is observed for the (la) and (Ic) data (see Table |3.3.2|) , however for (lb) 
(12^, p.b.c) there is instead a shift towards smaller sizes (see Table 3.3.2 and the increase 
of the number or instantons of size ~ 0.4 fm indicated on Fig. 13 by an arrow). We 
believe this effect to be mostly due to the instability of isolated instantons with p.b.c. (see 
section |2.4| ). To check this conjecture we left out of the ensembles (la,b,c) conhgurations 
with action S = Stt^. In Fig. |T^ we compare the original distributions (la,b,c) with the 


ones generated from them by removing S = Svr^ conhgurations (keeping the normalization 
of the original distributions), which we denote by (la’,b’,c’) respectively. We see that the 
boundary conditions dependence is reduced, supporting the conclusion that for L = 1.44fm 
this is the dominant hnite size effect. 

We can now go back to the discussion in section |2.4| about the strong hnite size ehects on 
the data ( 2 ) taken at L = 1.02 fm. As we have just discussed, physically relevant instantons 
have a typical size of 0.43 fm. The bulk of the size distribution is thus concentrated in the 
region where we have argued the ehect of the boundary begins to be strong and distort 
the instantons. L = 1.44 fm however seems to be just safe enough. A general conclusion 
we extract from this analysis is that to be free of hnite size ehects the lattice physical size 
should be at least 3 times the typical physical instanton size. 

The cutoh at p ~ 2.3a of our method eliminates lattice artifacts by construction. 
Another approach would consist in measuring the size distribution of all objects, physical 
and unphysical, then to separate the ones from the others by a scaling study. This is the 
approach attempted in , using “over-improved” cooling with various over-improvement 
coefficients. In that case the size distribution will be modihed from Fig. 0: it will 
show in addition a dominant, divergent contribution of lattice artifacts at small p. While 
this approach has the merit of extracting as much information as possible, the physically 
relevant piece is much harder to separate. 


3.4 Further properties of the Instanton Ensemble 

3.4.1 Instanton - Anti-Instanton Pairs 

Besides instantons and anti-instantons of physical sizes, the Monte Carlo simulation with 
Wilson action produces very many narrow peaks of both positive and negative topological 
charges. The structure at the level of “dislocations” disappears very soon in the cooling 
(~ 5 sweeps) independently of the lattice spacing, while well separated instanton - anti- 
instanton pairs of physical size may survive cooling for a while, although eventually they 
have to disappear due to their interaction (pairs are not minima of the action). 
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Because of the unphysical enhancement of short range structure (see for an evalua¬ 
tion of this effect) it is difficult to hx precisely the moment when annihilation of physically 
relevant I-A pairs starts. Judging, however, from our data, it seems that after 20 cooling 
sweeps the UV noise at the scale of the cut-off is strongly reduced, independently of {3. 
However it is not easy to design a scaling criterion for the I-A structure since it changes 
during cooling and moreover the speed of the changes depends on the lattice spacing and 
on the particular cooling algorithm. Our attitude has been to compare all quantities for 
different numbers of sweeps: whenever they are stable with cooling, as happens for instance 
with the charge and size distributions in Figs. 0, O, we are confident in considering them 


among the physical attributes of the original Monte Carlo configuration. Since above 20 
sweeps all the changes induced by cooling are mostly due to I-A annihilation, some infor¬ 
mation about the process of annihilation can be extracted by looking at the evolution with 
cooling of several quantities, like the total number of I’s and A’s or the smallest separation 
between instanton and anti-instanton. For this argument the stability of instantons under 
prolonged cooling, as ensured by the 5Li action, is essential. 

The presence of many close I-A pairs in the process of annihilation imposes certain limits 
on the reliability of our programs for identifying instantons (a local structure deformed by 
the annihilation may show many close peaks which should not be identified as separate 
instantons). This is taken into account to some extent by the self-duality cut but in some 
cases, mostly for small numbers of cooling sweeps, this may fail. Based on the comparison 
between different determinations of the instanton size for the 12^ and 12^ x 36 lattices (see 
section |3.3|) we give more than 90% conhdence to our “peak program” after 20 cooling 
sweeps or more. Since the peak assignment after 5 sweeps has a lower precision we chose 
not to present results based on it for such short cooling. 

What short cooling means from the point of view of I-A annihilation depends both on 
the lattice spacing and on the cooling algorithm. The annihilation of a pair at a certain 
physical distance is slower for smaller a since the effect of cooling a particular link extends 
at most a distance 3\/2a from it for bLi cooling. On Table |l| we give the number of peaks 
at various cooling stages, compared with the average action and absolute value of the 
topological charge. A comparison between the 12^ (a = 0.12 fm) and the 24^ (a = 0.06 
fm) lattices shows that on the latter about two or three times more sweeps are needed 
to reduce the number of pairs to the same level as on the former. While dislocations 
annihilate after a number of cooling sweeps roughly independent of /?, physical I-A pairs 
annihilate after a number of sweeps which increases with /3. Before 20 cooling sweeps for 
the 12^^ lattices and 50 sweeps for the 24^ lattice the abundance of close, annihilating pairs 
makes the description of the topological structure uncertain. To compare similar situations 
for the two lattices we present the results based on the peak analysis for 20 (50) sweeps 
and above. Note that the topological charge saturates much earlier (around 5 sweeps, see 
also Table |3.2|a) so the discussion above does not affect at all the determination of the 


^ Note that, even in the continuum, the problem of identifying all topological objects is not well-defined: 
classifying overlapping objects of positive and negative topological charges as an I-A pair (to keep) or a 
trivial fluctuation (to smooth out) becomes arbitrary as their separation becomes comparable to their size 
and the dilute picture breaks down. Different degrees of smoothing - as achieved at different stages of 
cooling - select in some sense different distance and overlap scales for the identihcation of the topological 
structure of the configurations. In our “short cooling” analysis we try to extract as detailed topological 
information as possible while still preserving a reliable instanton identification. 
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topological susceptibility. 



12 "^ t.b.c 

(la) 



12 ^ p.b.c 

. (16) 




f3 = 2A 




f3 = 2A 




sw 

(A^peaks) 

(^) 

(IQI) 

(Q^) 

(A^peaks) 

(5) 

(IQI) 

(Q^) 

20 

10 .6(1) 

8.42(7) 

1.73(7) 

4.5(3) 

10 .6(2) 

8.3(1) 

1.7(1) 

4.5(5) 

50 

5.1(1) 

4.48(6) 

1.69(7) 

4.4(3) 

5.08(8) 

4.34(4) 

1.6(1) 

4.3(5) 

150 

2.45(7) 

2.20(4) 

1.65(8) 

4.2(3) 

2.26(6) 

2.15(5) 

1.6(1) 

4.0(4) 

300 

2 .02(8) 

1.69(7) 

1.60(7) 

4.0(3) 

1.7(1) 

1.6(1) 

1.5(1) 

3.9(4) 


12 ^ X 36 

p.b.c. (1 

c) 


24^ t.b.c 

(3) 




P = 2A 




p = 2.6 




sw 

(-Speaks) 

(^) 

(IQI) 

(Q^) 

(-Speaks) 

(^) 

(IQI) 

(Q^) 

20 

31.7(4) 

25.3(3) 

2.9(2) 

13(1) 

63.1(5) 

30.7(7) 

1.7(2) 

4.4(8) 

50 

15.4(3) 

13.2(2) 

2 .8(2) 

12(1) 

15.7(5) 

12 .8(6) 

1.5(1) 

4.3(7) 

150 

6.6(3) 

5.9(2) 

2.7(2) 

12(1) 

6.7(3) 

6.1(3) 

1.6(2) 

4.2(7) 

300 

4.5(2) 

4.0(2) 

2 .6(2) 

11 (1) 

3.9(2) 

3.8(2) 

1.6(2) 

4.0(6) 


Table ^ Cooling behavior of the average number of peaks (iVpeaks), {S), (|Q|) and {Q^)- 
The (Ic) data for (iVpeaks), (S), (|Q|) should be rescaled by a factor and for {Q^) 

by a factor Ng/Nt in order to be compared with the others. 

In asking which are the typical features (size, location) of the physically relevant pair 
structure we can, to a certain extent, look for an answer by comparing the distributions, 
say, after 20(50) and after 300 cooling sweeps. Although between the two cooling stages 
about 80% of the (anti-)instantons have annihilated (see Table |l]) the shape of the size 
distributions is very similar. This is clearly visible on the normalized distributions of 
Fig-H and suggests that the size distribution given there holds also for the instantons and 
anti-instantons appearing in pairs, which have not annihilated in the very early stages of 
cooling. Since this distribution scales correctly with the cut off this is simultaneously a 
consistency check that no unphysical structure was present at 20(50) sweeps. This bound 
is surely an overestimate, as has been argued above. We have, however, no reliable way of 
characterizing the physical structure bound in close pairs which annihilate between 5 and 
20(50) sweeps. 

3.4.2 Distance Distributions 

In Fig. 1^ we present distance distributions for alike and opposite charge objects. The 
solid line is the 4-dimensional volume factor, that is it corresponds to a homogeneous 
distribution of instantons over the lattice (the shape is due to the hnite size of the hyper- 
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cubic, periodic lattice). The figure shows little difference with a homogeneous distribution. 
After 300 sweeps most of the I-A pairs have disappeared and are thus not included in the 
figure. 

In Fig. 0 we present distributions of distances to the closest alike and opposite charge 
objects. The center of the distribution shifts toward larger distances with longer cooling as 
a result of pair annihilation. For the 24^ lattice we observe also a small peak in the alike 
distribution at short distances (~ 3a) and short cooling, which disappears with further 
cooling. A much stronger peak at small distances (about 1 to 2a) has been observed by 
|2l|| . In our opinion structure at the scale of the cut-off is unphysical and in fact the cuts 
we impose in the peak assignment (see section |3.1|) take care of most of it. The remnant 


at 3a may have physical significance. We believe however that it is to be attributed to 
a bad assignment of peaks. A wide, deformed instanton could also trigger such an effect, 
but it should not depend then on the amount of cooling. We think that the annihilation 
of I-A pairs is the dominant effect: in the course of pair annihilation, transient deformed 
structures can appear which lead to multiple alike-charge peaks. These then automatically 
produce a temporary increase in the small distance part of the alike distribution, an effect 
which is removed by further cooling. The cuts we imposed from self-duality take care of 
most of this effect but still part of it has remained. 

An interesting information provided by Table is that the number of pairs, Apairs, 
observed after short cooling - 20 sweeps for the 12^ lattice and 50 sweeps for the 24^ lattice 
- is significantly larger than the value corresponding to a binomial distribution with equally 
probable, uncorrelated I and A events, which would give 


2{N^„) + (IQI) = {Afp..ta> = (Q">. 


(23) 


Simultaneously we observe (*F)/(Apeaks) < 1 (see also Fig. in next section). The distri¬ 
butions at 20(50) sweeps seem in fact compatible with binomial laws for uncorrelated I, A 
(with equal probabilities such that {Nj + Na) = (Q^), see eq. (^3|)) and I-A events (ac¬ 
counting for the surplus of I-A pairs, (AApairs) = ((Apeaks) — (Q^))/2), see the histograms 
in Fig. 00. Due to the annihilation of pairs during cooling this distribution is modified 
and after long cooling evolves towards a binomial distribution for only one type of charges 
(as after 300 sweeps for lattices (la,b)). These facts taken together are compatible with an 
I-A interaction, which would be responsible both for enhanced pair creation in the Monte 
Carlo production and for pair annihilation during cooling. This interaction is surely de¬ 
pendent on the overlap between opposite charge objects. We can characterize the overlap 
for I-I or I-A pairs as 

^ ^ Pi + P2 

_ ri2 

^ For comparison with the actual MC data, we generated spatial configurations of charges with the 
same susceptibility (Q^) and number of peaks (A^peaks) as follows. Pairs of charges are placed on the 
lattice with an average density dp, at locations x and y chosen independently with uniform probability. 
They do not contribute to the susceptibility, but give rise to 2dpV peaks on average. In addition, single 
charges are added with average density ds, at random locations, contributing an additional dsV to the 
number of peaks, so that (A^peaks) = {ds + 2dp)V. If the sign of these single charges is chosen randomly, 
the susceptibility will be (Q^) = dgV. This is appropriate if (Apeaks) > {Q^), ie. if there is a surplus of 
pairs over a binomial distribution. If instead there is a defect of pairs, as for long cooling, the sign of the 
single charges is chosen “en bloc” for each configuration. The susceptibility is then (Q^) = dgV -I- (dsV)'^. 
In both cases, the location of each charge is uncorrelated with that of the others. 
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where pi and p 2 are the sizes of the objects in the pair and ri 2 their mntnal distance. We 
obtain for the average over (la,b) lattices after 20 sweeps (O) ~ 1.7 for I-A pairs and 
(O) 1.8 for I-I pairs. On the (3) lattice after 50 sweeps (O) ~ 1.7 for I-A pairs and 
(O) 2.0 for I-I pairs. The overlap we observe is rather large, suggestive of a percolating 

polymer structure rather than a dilute gas. 

Since the action is lowered by the interaction, which increases with the overlap, we 
would expect increased production of opposite pairs with larger overlap. Surprisingly 
however, there is no indication after 20(50) cooling sweeps of a larger overlap for opposite 
charge than for alike charge pairs, which we would naively expect to result from enhanced 
production of close I-A pairs in the original Monte Carlo data. Hence we can only say that 
the I-A interaction favors conhgurations with a surplus of I-A pairs over a naive binomial 
distribution, but with a rather homogeneous, uncorrelated distribution of I’s and A’s over 
the lattice. Since cooling also preferentially destroys pairs with large overlap we cannot 
exclude, however, that close I-A pairs were more abundant in the uncooled conhgurations. 


3.4.3 Fractional Charge Instantons 


It has been recently pointed out that instantons with topological charge 1/2 might 
play a relevant role in the SU(2) vacuum. The presence of such structure could modify 
the general picture of the instanton ensemble. In particular it will change the statistical 
relation between susceptibility and number of peaks, eq. (|^), since such objects carry only 
half the topological charge and action of ordinary instantons. Part of the analysis in 


to detect them was based on a determination of the number of peaks in the energy density 
compared with the average value of the action. If instantons of charge 1 are the dominant 
structures in the cooled Monte Carlo conhgurations one expects that the average number 
of peaks will coincide with (S). Instead in Ref. 
found. 
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an excess in the number of peaks was 
We have pursued a similar analysis over our conhgurations. In Table we have 


given the number of peaks compared with {S) as a function of cooling sweeps. There is an 
excess of Apeaks over {§) but what we observe amounts in general to 20% or less, at short 
cooling, and decreases fast with cooling, in such a way that for 300 sweeps both numbers are 
roughly compatible within errors. This indicates that if such structure was there initially 
it should have appeared mostly in opposite charge pairs which have annihilated after long 
cooling. It is particularly difficult to distinguish a Q = ±1/2 pair from an ordinary I-A pair 
and the above ehects could be due also to the latter. The presence of several topological 
objects in our physical box prevents the integration of the topological charge density over 
all space: when these objects are large compared to their separation and form a dense gas, 
the total charge carried by each object becomes ambiguous. After 300 sweeps we looked 
directly into some of the few conhgurations for which there was an excess in the number of 
peaks. We found some suggestion of the presence of Q = 1/2 objects, but we do not know 
a good way of identifying and separating them further from some other kind of structure 
like wide instantons, a problem already mentioned in 
discussion are rather close to the uncertainty level of 
identifying peaks. 

In Fig. |T^ we present the distribution of /ip = *§/Apeaks- If the Q = 1/2 were the 
dominant objects the points would concentrate mostly around /ip = 1/2 while for ordinary 
instantons they will concentrate around /ip = 1. Our results show a clear tendency towards 


In any case, the effects under 
10 % we assign to our programs for 
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/ip — 1. 

On the basis of these results we conclude no clear indication for the presence of such 
objects in our MC conhgurations. Note, however, that due to the periodicity of our bound¬ 
ary conditions in spatial directions, fractional charge instantons cannot appear alone. It 
is possible that to see a clear effect of these objects under these conditions larger volumes 
than the ones used in our simulations are needed (the ones employed in are almost 
twice as large as ours). This is an interesting proposal which deserves further investigation. 


3.5 Density - Density Correlations 

Using the improved charge operator and action allows to measure meaningful density - 
density correlations even after a small number of cooling sweeps. The extraction of the 
instanton features from these data is more involved and therefore matter for a separate 
investigation. We show here part of these data on Figs. and ^ only to indicate the com¬ 
patibility with results of the previous analysis. We note by q+{x) and q-{x) max(g(x),0) 
and min(g(a;),0) respectively. 

Several features are noticeable: 


• For short cooling (g(O)g(r)) goes smoothly to zero as r increases, as expected, in¬ 
dicating a rather homogeneous distribution of instantons and anti-instantons. For 
long cooling the tail represents a more or less homogeneous distribution of a few 
instantons or anti-instantons, since practically all pairs have meanwhile disappeared. 


• As r —>■ cx), lim(g(0)_|_g(r)+) = lim(g(0)_g(r)_). This is a consequence of the evenness 
of the ^-distribution. For short cooling and as a consequence of the homogeneous 
distribution of instantons and anti-instantons we see in addition hm(g(0)+g(r)+) = 
—lim(g(0)+g(r)_) = ilim(|g(0)||g(r)|), within statistical errors. For long cooling, 
however, as a consequence of the disappearance of most I-A pairs hm(g(0)+g(r)_) 

0, and lim(g(0)+g(r)+) = ilim(|g(0)||g(r)|). 


The peak of the q — q correlations near the origin is caused by the typical size of the 
instanton. It can be used to extract information about the instanton size |^. For 
comparison we present (dotted line in Fig. 13 ) the distribution obtained with one 
instanton obeying the continuum ansatz eq. (|]), with = pm = 0.43 fm, rescaled to 
match the measured value of the {q{0)q{r)) correlation at the origin. The moderate 
agreement shows the limitations of a one-instanton ansatz in interpreting the density 
correlations (see below for a more realistic ansatz). 


• Comparing the q — q and |g| — |g| correlations we see clearly the minimal distance 
between opposite charges, in agreement with vanishing of the g+ — g_ correlation 
below this distance. 


To help us interpret our Monte Carlo data on density-density correlations, we have 
constructed some “synthetic” data from the simplest ansatz: we place on the lattice in¬ 
stantons obeying the continuum ansatz (with periodicity effects accounted for by adding 
the 80 nearest replicas). These “synthetic” instantons do not interact. They are distributed 
homogeneously (with no spatial correlations) as described in section |3.4.2| , with densities of 
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instantons and of I-A pairs adjusted to match the susceptibility and the number of peaks 
measured in the Monte-Carlo. The size distribution for these continuum instantons is cho¬ 
sen according to the measured size distribution. We can then compare our density-density 
correlations with those produced by this parameter-free ansatz. The ansatz results are 
given by the solid curves in Figs. |^. The agreement with the actual Monte Carlo data 
is excellent (it becomes almost perfect if one allows for an overall normalization factor, 
to account for, e.g., the small mismatch of (Q^) and (iVpeaks))- This validates a posteriori 
several aspects of our approach: deviations of our lattice instanton fluid from a linear super¬ 
position of continuum ansatz are very small; instanton sizes have been determined reliably. 
Only the correlations of the action density at short cooling are not so well reproduced by 
our non-interacting ansatz. This could be expected, since our ansatz assigns an action of 
per charge, whereas the action we actually measure at 20 cooling sweeps is smaller 
(see Table 3.4.1). This is our hrst signal for an instanton-anti-instanton interaction. Its 
influence on the spatial charge distribution, however, is too weak for us to detect directly 
here. 

The features of the density-density correlations listed above are preserved when going 
to smaller scales, as indicated by comparing Figs. 0 and 1^ (notice that the normalizations 
are different, since the number of peaks is larger on the 24^ lattice than on the 12^ lattice, 
see Table 1^. It would be very interesting to compare these correlations with the predictions 
of phenomenological models, and to observe the effect of temperature and of dynamical 
fermions. We are starting work in these directions. 

The correlation {q{0)q{r)) of the topological charge density cannot be directly compared 
with the density-density correlation of a classical fluid made of positively and negatively 
charged particles. We do not deal with point-like or hard-core particles and the q — q 
correlation smears the information about the center of the instanton. We could apply a 
de-convolution, to reconstruct the correlation for an instanton to be at r, given the center 
of another one at the origin. Instead we measure such correlations directly, by placing a 
charge Q = ±1 at the location of the peak of q. The correlations {q{0)q{r)) thus obtained 
after this “lumping” of the topological charge are displayed in Fig. |^. Since there are few 
peaks at small distances (in fact we identify peaks at distances up to 1 in each direction), 
there is no contribution to these correlations here. By comparing these data with those 
for the full densities we clearly see the effect of the size of the instanton at short distances. 
The q — q correlation is zero both at small and at large distances, corroborating the smooth 
decrease {q{0)q{r)) with r, and conhrming the overall absence of structure in our instanton 
fluid. The small difference between the — g_ and the — g_|_ in the intermediate region 
may be a signal of pair annihilation. This would explain also the slight increase of |g| — |g| 
in the same region. For further conclusions the data are too noisy yet. 

As a hnal check for the appearance of regular (“crystalline”) structures, we analyzed 
the Fourier Transform of the density-density correlations. The spectrum does not show the 
typical peaks produced by a crystal. In agreement with our observations from the distance 
distribution, we conclude that there is a strong indication for a homogeneous, spatially 
uncorrelated distribution of instantons over the lattice. 
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4 Conclusions 


The goal of our study is to identify classical objects (instantons) on the lattice. For this pur¬ 
pose, it is appropriate to use a classically improved lattice action and topological charge 
operator. Our 5-loop operator has a tree level discretization error starting with 0{a^), 
which is hne-tuned to minimize violations of scale-invariance. We use this operator both 
for measurements, and for cooling, which rapidly isolates the topological components of 
the QCD vacuum. With this improved cooling algorithm we have thus a straightforward 
method to remove UV noise and dislocations without modifying the topological structure 
at physical scales (in so far as it is not affected by hnite size effects or pair annihila¬ 
tion). The excellent scale invariance properties of our cooling action allows us to identify 
unambiguously the physical structure of the original conhguration which carries the net 
topological charge, since instantons are completely stable under any amount of cooling. 
A scaling analysis of the results at various levels of cooling shows that already before 20 
improved cooling sweeps UV noise and dislocations have disappeared and only physical 
structure remains, including I-A pairs. The improved cooling, being a minimization pro¬ 
cedure, does not preserve I-A pairs, but since there are no other modihcations of physical 
instantons during cooling we can estimate the contribution of the pairs by comparing re¬ 
sults at different cooling stages. The improved charge operator gives smooth density data 
already on still rough conhgurations, data which can be used also for correlation analysis. 
The total charge comes out an integer within (9(1%) after 5-10 sweeps and remains stable 
thereafter. 

The dislocation cut off is hxed in lattice units (po — 2.3a) and therefore will not affect 
the physical scales for a small enough (such that po falls below these scales). The hnite 
size effects begin to affect instantons of sizes above L/3. Therefore we should take care to 
use lattices larger than three times the typical instanton size. 

The identihcation and description of the structure in the charge and energy densities 
can successfully be based on the continuum ansatz for instantons, especially if one takes 
into account anisotropy effects induced by quantum huctuations and clustering effects. 
Uncertainties in the identihcation procedure ahect either large instantons deformed by 
hnite size ehects or unstable structures like I-A pairs in late stages of annihilation. 

Using improved cooling we have analyzed the topological content of the SU{2) theory, 
especially the topological susceptibility and the instanton size and distance distributions. 
Since the decisive criterion for the physical signihcance of the results is the correct scaling 
with the cut oh in approaching the continuum limit, we worked at diherent scales. We 
checked that under a rescaling of a by a factor of two the susceptibility and the relevant 
part of the size distribution both scale very well. We also verihed that the lattices we 
used (size 1.4 fm, 12"^ and 24^) were large enough not to introduce hnite size ehects for the 
relevant physical sizes. Our results for the susceptibility and the size distribution are stable 
under cooling, which only ahects the distance distribution as the result of pair annihilation. 

The value of the susceptibility, 200(15)MeV, agrees well with the Witten-Veneziano 
formula. The most likely instanton size turns out to be 0.43(5) fm, with a full width of about 
0.47 fm. This size is larger than the one proposed for the instanton liquid models p ~ l/3fm 
0. The distance distribution is similar to that of randomly placed objects, and there is no 
indication for a crystalline structure. The structure of conhgurations at scales above the 
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threshold po ~ 2.3a - which we argued to be the physically relevant ones in what concern 
topological excitations - appears to be a random mixture of instantons and anti-instantons 
of sizes compatible with an effective semi-classical distribution and with a homogeneous 
spatial distribution. With short cooling (20 (50) sweeps for the 12'^ (24*^) lattices), however, 
we observe more pairs than would be compatible with uncorrelated binomial distributions 
for instantons and anti-instantons, suggesting that due to pair interaction the production 
of opposite charge pairs is enhanced. No significant correlation between the locations of A 
and I objects is seen here, but we cannot exclude the presence in the original configurations 
of closely bound I-A pairs which annihilated early during the cooling. The relative overlap 
eq. (1^) of about 1.5-2 among closest objects makes the dilute description questionable 
at the scale of short cooling - 20 sweeps on the (la,b) lattices. The density of instantons 
there is about 2 — 3/fm^ and the fraction of volume they occupy, counting \p\,[ for each, 
is about 30 — 50%. The long cooling regime is characterized by the absence of pairs and 
by a binomial distribution for the surviving type of charges, with a density set by the 
scaling susceptibility. Note that these numbers depend on the criterion for identifying 
instantons in the presence of I-A pairs, which is subject to an ambiguity existing already 
in the continuum. Finally, there is no clear indication for fractionally charged instantons, 
but for this analysis larger volumes might be necessary. 
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sweeps ^ sweeps 

Figure 1: Evolution of: (a) action, S, and (b) size, Ppeak (0))! with the number of 

cooling sweeps for various choices of the improvement coefficient C 5 (i.e. C 5 = 1 / 10 , C 5 = 0 
and C 5 = 1/20 for 3Lz, ALi and 5Li respectively). 



Figure 2: Evolution of instanton size under 5Li improved cooling for various initial con- 
hgurations. Notice the stability of the size with cooling for instantons above the threshold 
Po = 2.3a. 
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Figure 3: Dependence of the 5Li action (squares) and Wilson action (circles) on the 
instanton size. The triangles are data coming from Monte Carlo conhgurations generated 
with the Wilson action and cooled with 300 sweeps of 5Li cooling. The vertical dashed 
line indicates the location of the threshold po = 2.3. The dotted curve is the value of the 
Wilson action for an instanton of size p calculated at tree level to it shows the 

deviation from scale invariance of the Wilson action. 
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Fignre 4: Behavior nnder cooling of conhgurations with one instanton. We show the 
dependence on lattice geometry and boundary conditions. Sqnares and circles correspond 
respectively to 5Li improved cooling and over-improved cooling with e = — 1. We also 
present for comparison a 2-instanton conhguration. 
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Figure 5: Values of the topological susceptibility for the various lattices: t.b.c. 
(la)(squares), (2)(stars) and (3) (crosses), and p.b.c. (lb) (triangles) and (Ic) (circles) 
- see Table ^T| . In lattices (la,b,c) we have suppressed the error bars to make the £g- 
ure clearer, the bars are an indication of the magnitude of the error for, from right to left, 
(la,b,c). In the upper plot we present the dependence of the topological susceptibility with 
lattice spacing after 20 and 300 cooling sweeps. The point at a = 0.12fm is an average 
over lattices (la,b,c). 
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Figure 6: Charge distribution after 20 and 300 bLi cooling sweeps for the lattices: (la) 
(squares) and (lb) (triangles) - see Table 1^ 
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Figure 7: Upper plot: size distributions after 20 cooling sweeps extracted from peak 
size ppeak, eq. (0). and “center and nearest-neighbor” fit pc, eq. ([T^). We also show 
the histogram of excentricity Cc, eq. ( 0 ) and (|T^, after 20 and 300 cooling sweeps. 
The continnous, dotted and dashed histograms correspond to lattices (la), (lb) and (Ic) 
respectively (see Table ^T|) . Lower plot: dependence of the excentricity Cc on instanton 
size, lattice geometry and boundary conditions. 
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Figure 8: Comparison between different size definitions for: upper plot: 20 cooling sweeps; 
lower plot: 300 cooling sweeps. Data are taken on (la,b,c) lattices. 
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Figure 9: Comparison between instanton size distributions extracted from Ppeak and p\. 
Data are obtained by cooling for 300 sweeps configurations obtained from a Monte Carlo 
simulation on a 12^ lattice, with {3 = 2.5 and twisted boundary conditions. N{p) denotes 
the number of instantons with size p. 



Figure 10: Histogram of normalizations (see eq. (||)). Upper plot: 20 cooling sweeps; lower 
plot: 300 cooling sweeps. 
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Figure 11: Size distributions. Crosses, squares, triangles, circles correspond respectively 
to lattices (3), (la), (lb), (Ic) (see Table |0|) , the dotted vertical line is the small size 
threshold po = 2.3a. The curve in the hgure corresponds to the £t eq. (PH) . 
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Figure 12: Size distributions before and after subtracting configurations with ^ = 1. 
Squares, triangles and circles correspond to (la), (lb) and (Ic) respectively (see Table 
0 ), the dotted vertical line is the small size threshold po = 2.3a. 
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Figure 13: Distribution of separation between instantons. Circles and squares represent 
separation between alike and dislike-charged peaks respectively. The solid-line histogram 
is the 4-dimensional volume factor corresponding to a homogeneous distribution of (anti- 
)instantons over the lattice. The upper plot is an average over (la) and (lb) lattices. The 
lower plot corresponds to the (3) lattice. 
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Figure 14: Distribution of separation to the nearest alike-charged (circles) or dislike-charged 
(squares) peak. The upper plot is an average over (la) and (lb) lattices. The lower plot 
corresponds to the (3) lattice. The histograms in the hgure represent simulations from 
homogeneous binomial distributions as described in section Eg. 
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Figure 15: Distribution of /ip = S'/iVpeaks with iVpeak, the number of peaks in the energy 
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Figure 16: Density-density correlations after 20 cooling sweeps for lattice (lb). We plot 
correlations of the dimensionless densities obtained from the continuum ones by multiplying 
with the physical volume of the lattice. In the upper part of the plot we show (s(O)s(r)), 
(g(O)g(r)) and (|<?(0)| |g(r)|) denoted as ss, qq and |g||g| respectively. In the lower part 
of the plot we have magnihed the scale to show a comparison between {q{0)±q±{r)) and 
(g(0)+q'_(r)) {q±q± and q+q- respectively in the hgure). The solid curves have been ob¬ 
tained from a homogeneous binomial distribution of instantons obeying the continuum 
ansatz as described in sections |3.4.2|, |3.4.1. 
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Figure 17: Density-density correlations after 300 cooling sweeps for lattice (lb). The dotted 
line represents the correlation obtained from one instanton of size 0.43fm rescaled to match 
the measured value of the {q{0)q{r)) correlation at the origin. The solid curves have been 
obtained from a homogeneous binomial distribution of instantons obeying the continuum 
ansatz as described in sections |3.4.2| , |3.4.1| (on the upper plot only the curve for qq appears 
since due to the small number of I-A pairs the ss and \q\\q\ curves practically coincide with 
this one). Notation as in Fig. 
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Figure 18: Density-density correlations at 50 cooling sweeps for 20 configurations from 
lattice (3). The solid curves have been obtained from a homogeneous binomial distribution 
of instantons obeying the continuum ansatz as described in sections |3.4.2 , 
as in Fig. [TO. 


3.4.1. Notation 
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Figure 19: Density-density correlations after 20 cooling sweeps for lattice (lb) obtained by 
placing a charge Q = ±1 at the location of the peaks of q{x) and zero everywhere else. 
Notation as in Fig. 


16. 
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Figure 14: Distribution of separation to the nearest alike-charged (circles) or dislike-charged 
(squares) peak. The upper plot is an average over (la) and (lb) lattices. The lower plot 
corresponds to the (3) lattice. 














